TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

电子设备吊舱瞬态热载荷分析与计算

瞬态热载荷的分析与计算是合理设计电子设备吊舱环境控制系统的基础.本文建立了吊舱动态热平衡方程与壁面导热方程,并使用有限差分方程求解吊舱瞬态热载荷.分析了飞行状态及吊舱热容量对吊舱瞬态温度的影响,提出了"蓄冷"节能的设计概念并据此合理地确定了吊舱环境控制系统的设计热载荷.     

弧形反射罩几何形状与到达热流密度均匀性关系研究

文章在推导出兰贝特定律的平面形式的基础上,利用Monte Carlo方法分析了弧形红外灯反射罩的各项参数及灯间距对辐照均匀性的影响.     

飞机风挡鸟撞动响应分析方法研究

建立了鸟体模型和风挡玻璃破坏准则,形成了一套完整的风挡鸟撞动响应分析方法。并对某型号风挡进行了鸟撞动响应分析,得出了鸟的运动轨迹和风挡的位移、应力、应变响应以及风挡玻璃的鸟撞临界速度,经与试验结果比较,两者吻合较好。     

红外加热笼边缘效应对卫星表面热流均匀性的影响研究

文章研究了红外加热笼边缘效应对卫星表面到达热流密度均匀性的影响,主要通过增加边缘加热笼并将其翻转一定角度的方法,对红外加热笼的设计进行优化,以改进卫星表面的到达热流密度均匀性.建立了变角度情况下蒙特卡罗仿真计算的几何模型,对3种情况进行了计算和分析,得到了对红外加热笼的边缘效应进行优化设计的方法,可以用于指导多种情况下的红外加热笼设计.     

航空发动机高压涡轮叶片表面红外发射率测量及应用

在使用RotamapⅡ红外测试系统前需要已知被测涡轮叶片的表面发射率。阐述了1种使用电涡流加热涡轮叶片,从而获得较大温度范围内的涡轮叶片表面发射率的方法。试验测试了某型航空发动机的新、旧2种高压涡轮叶片的表面发射率。试验结果表明:不同温度下叶片的发射率数值会发生变化,新叶片的表面发射率随温度的升高而明显增大,当叶片表面发生严重氧化后发射率变化较小,在±0.02之间。采用RotamapⅡ系统测试某型发动机高压涡轮叶片温场时,发射率数值可以取0.893,其测温误差小于被测物体温度的1%。同时结合试验得出的发射率造成的测温误差曲线可以对红外测温结果进行修正。     

民用飞机风挡表面对流换热系数的校核

适航条款要求在结冰条件下,应确保风挡表面不能结冰,防止影响飞行员视界。因此应确定结冰条件下具有足够的电加热功率使风挡表面温度高于冰点(零度),而影响加热功率的关键因素是外表面的对流换热系数。为了获得准确的对流换热系数,在低速风洞中测量了不同工况下风挡表面对流换热系数,并利用CFD方法建立了一个仿真风洞模型,计算试验工况下的风挡表面对流换热系数,通过测量值来校核CFD的计算值,获得一个修正方法,最终使用这个修正方法计算获得其他工况下的对流换热系数。     

驾驶舱透明件安装形式分析

主要描述了世界上民用飞机驾驶舱透明件(主要是风挡)所采用的三种安装形式,对这三种安装形式的优缺点进行了多方面的分析和对比,并从结构设计的角度进行了客观评价。     

基于再入轨迹和气动热环境的返回舱烧蚀研究简

 针对再入全过程合理预测热防护罩表面材料烧蚀深度和温度的动态变化问题,提出融合再入轨迹、气动热以及Newton-Raphson和三对角矩阵算法（TDMA）构建动态烧蚀的方法。该方法建立直入式和跳跃式三自由度再入轨迹,应用修正的牛顿流体理论估算气动参数,以及修正的Fay-Riddell和Sutton-Grave理论计算驻点区域的热流密度,利用一维非线性热传导方程模拟了热防护材料的烧蚀过程。仿真结果表明：此方法实现了再入全过程热防护材料烧蚀深度和温度连续动态变化的预测,同样适用于更为复杂结构飞行器的动态烧蚀预测,与热平衡积分法（HBI）相比其结果可靠合理,为进一步优化热防护系统（TPS）提供了一定的参考依据。     

Numerical investigation of windage heating within shrouded rotor-stator cavity system with central inflow

The rotating disk surface temperature rise due to windage heating effect by numerically modeling the turbulent flow within a rotor-stator cavity which is available with a peripheral shroud and imposed through airflow was dealt with. The windage heating may be defined as viscous friction heating caused by relative velocity differences across the boundary layers between the fluid and the rotating disk surface. The kinetic energy dissipation process could transform the rotating shaft power into thermal heating. Commercial finite volume based solver, ANSYS/CFX was employed to numerically simulate this physical process by using the shear stress transport (SST) turbulence model. CFD results include the rotating disk surface temperature axial distribution and tangential velocity distribution of the fluid domain. The velocity difference between the result obtained by particle image velocimetry (PIV) experiments and CFD simulation are within 5%. The adiabatic disk temperature rise can be calculated by the tangential velocity of disk and fluid in large gap ratio and turbulent parameter. CFD temperature distribution results and those estimated via velocity differences are within 10%.